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The entanglement entropy (von Neumann entropy) has been used to characterize the complexity 
of many-body ground states in strongly correlated systems. In this paper, we try to establish a 
connection between the lower bound of the von Neumann entropy and the Berry phase defined 
for quantum ground states. As an example, a family of translational invariant lattice free fermion 
systems with two bands separated by a finite gap is investigated. We argue that, for one dimensional 
(ID) cases, when the Berry phase (Zak's phase) of the occupied band is equal to n x (odd integer) 
and when the ground state respects a discrete unitary particle-hole symmetry (chiral symmetry), 
the entanglement entropy in the thermodynamic limit is at least larger than In 2 (per boundary) , 
i.e., the entanglement entropy that corresponds to a maximally entangled pair of two qubits. We 
also discuss this lower bound is related to vanishing of the expectation value of a certain non-local 
operator which creates a kink in ID systems. 



I. INTRODUCTION 

One of the most distinctive features of quantum phases 
of matter is that they are not completely characterized by 
their pattern of symmetry breaking (order parameters of 
some kind), which is in a sharp contrast to classical sta- 
tistical systems. Instead, quantum ground states should 
be described by their pattern of entanglement such as 
topological or quantum order. Q] However, beyond some 
simple textbook examples, e.g., a system of two coupled 
5 = 1/2 spins (qubits), we do not have many intuitions 
about quantum entanglement hidden in many-body wave 
functions. In a recent couple of years, the entropy of en- 
tanglement (von Neumann entropy) [2] 

S A = -tT A p A hip A , p A =tr B \y)m, (1.1) 

has been used to measure how closely entangled (or how 
"quantum") a given ground state wave function \ty) is. 
Here, the total system is divided into two subsystems 
A and B and pa is the reduced density matrix for the 
subsystem A obtained by taking a partial trace over the 
subsystem B of the total density matrix p = \^)(^\. 
This quantity is zero for classical product states whereas 
it takes a non-trivial value for valcncc-bond solid states 
(VBS), or resonating valence bond states (RVB) of quan- 
tum spin systems, say. Recently, the entanglement en- 
tropy at and close to quantum phase transitions in low- 
dimensional strongly correlated systems has been used as 
a new tool to investigate the nature of quantum critical- 
ly- HH SIS II Even though one can tell different 
quantum phases from the scaling of the entanglement 
entropy, it is still not completely understood what kind 
of information we can distill from the von Neumann en- 
tropy, other than that contained in conventional correla- 
tion functions. 

On the other hand, a phase degree of freedom is also 
a specific feature of quantum mechanics. Indeed, Berry 
phases |lOj associated with (many-body) wave functions 



in solid states are related to several interesting quan- 
tum phenomena which have no classical analogue. Prob- 
ably, it is best epitomized by the Thouless-Kohmoto- 
Nightingale-Nijs (TKNN) formula in the integer quan- 
tum Hall effect (IQHE) [Hill, 

in which gapped quan- 
tum phases are distinguished by an integral topological 
invariant originating from winding of the phase of wave 
functions. In addition to the IQHE, the Berry phase 
also appears in the King-Smith- Vandcrbilt (KSV) for- 
mula [13 . Il4j of the theory of macroscopic pola rization, 
and its incarnation in quantum spin chains |l5lll6l |. and 
so on. An observable consequence of the non-trivial Berry 
phase is the existence of localized states at the boundaries 
when we terminate a system with boundaries. |l7lll8lll9l | 

It is then tempting to ask if there is, if any, a connec- 
tion between these two paradigms in quantum physics, 
namely, entanglement and the Berry phase. In this pa- 
per, we discuss this issue by taking a family of transla- 
tional invariant lattice free fermion systems in d dimen- 
sions as an example. We bipartition the system into two 
subsystems A and B by introducing (e? — l)-dimensional 
flat interfaces. Within this setup, we can reduce the cal- 
culation of the entropy to that in a one-dimensional sys- 
tem by the (d — l)-dimensional Fourier transform along 
the interface. We assume the existence of a finite energy 
gap to above ground states which is inversely propor- 
tional to the correlation length, m ~ (,corr (when mea- 
sured in the unit of the band width). Furthermore, for 
simplicity, we consider the case in which there are only 
two bands that are separated by a gap. 

In this paper, we consider the Berry phase associated 
with a response of a quantum ground state to a contin- 
uous twist of the boundary condition. For the case of 
free lattice fermion systems, for which a ground state is 
given by a filled Fermi-Dirac sea, this Berry phase is a 
phase acquired by an adiabatic transport of the Bloch 
wave functions in the momentum space and also called 
Zak's phase. (2^ Physically, it is related to macroscopic 
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polarization of the Fermi-Dirac sea. |ia A beauty of 
the simple two-band example that we discuss is that the 
Berry phase for the quantum ground state can be easily 
computed and visualized, following the pioneering work 
by Berry (See SecHTJand Fig. [JJ below.) 

With these setups, wc will demonstrate that taking the 
partial trace over a subsystem corresponds to creating 
boundaries in a system. Two contributions to the en- 
tanglement entropy will be then identified. The first one 
is of type already discussed in a flurry of recent works 
focusing on detection of quantum critical points. This 
contribution to the entanglement entropy is largely con- 
trolled by the correlation length £ corr . For example, in 
one-dimensional (ID) many-body systems close to criti- 
cality the entanglement entropy obeys a logarithmic law 
Sa ~ A(c/6) ln£ corr /a where c is the central charge of 
the conformal field theory that governs the criticality, a 
the lattice constant, and A is the the number of boundary 
points of A. HQ 

On the other hand, the second contribution to the en- 
tropy comes from the localized boundary states of the 
correlation matrix that exist when the Berry phase of the 
ground state wave function is non- vanishing. Especially, 
when the Berry phase is equal to 7r x (odd integer) and 
when the ground state respects discrete symmetries of 
some sort, the localized boundary states are topologically 
protected as discussed in Refs. For this case, we 

will show that the contribution from the boundary states 
to the von Neumann entropy is In 2 per boundary, i.e., the 
same amount of entropy carried by maximally entangled 
pair of two qubits. We will also illustrate, by taking a 
specific limit, that when 7^0, the von Neumann entropy 
from the boundary states is that of partially entangled 
qubits. 

We also discuss that the In 2 contribution to the von 
Neumann entropy is related to vanishing of the expecta- 
tion value of a certain non-local operator which creates 
a kink in ID systems. This connection between the en- 
tanglement entropy and the kink operator is, in flavor, 
similar to discussions in Refs. Q, 0| in which the entan- 
glement entropy is expressed as the expectation values of 
twist operators in conformal field theories. 

The rest of the paper is organized as follows. In Sec. ITU 
we start our discussions with ID translational invariant 
Hamiltonians with two bands separated by a finite gap. 
The Berry phase is introduced as an expectation value 
of a specific non-local operator that twists the phase of 
wavefunctions. We then discuss its connection to the 
entanglement entropy by making use of the correlation 
matrix. The calculation of the entanglement entropy is, 
in general, a rather difficult task at least analytically. 
Furthermore, the Berry phase contribution to the en- 
tropy might not be of perturbative nature. We thus con- 
sider two limiting situations. In Subscc. lll Bl we take the 
limit of the small correlation length £ corr <C 1 and zero 
band width. In this specific limit, we can express the 
entanglement entropy as a function of the Berry phase 7. 
We next focus on cases with a discrete unitary particle- 



hole symmetry (chiral symmetry) in Subsec. Ill CI Ex- 
cept for requiring the chiral symmetry, any parameters of 
the Hamiltonian (the band structure) can be arbitrary. 
Once we impose the chiral symmetry, the Berry phase 
7 can take only discrete values, integral multiple of k. 
We then show when 7 = it x (odd integer), the entan- 
glement entropy is bounded below as Sa > 2 In 2. In 
Sec. IIIII we relate the lower bound of the entropy at 
7 = 7T x (odd integer) to the vanishing of a non-local op- 
erator that creates a kink. In Sec. II VI these discussions 
are applied to a higher dimensional example, a 2D super- 
conductor with non-zero TKNN integer. We conclude in 
Sec. El 

II. ID TWO-BAND SYSTEMS 

We start from the following ID translational invariant 
Hamiltonians with two bands separated by a finite gap, 

PBC / A \ 

H = E ■■<■-'-■ H *-*' = { A> i _) , ' ^ 

Here, a pair of fcrmion annihilation operators = 
(c + ,c_) x is assigned for each cite, x,x' = 1, ■•■ ,N, 
and the hermiticity of TL implies t L ^ x - x i = t* x ,_ x and 

Ax-x' = {^x'~x)* f° r t = ±. We impose the periodic 
boundary condition (PBC) on the ID lattice. In spite of 
its simplicity, this Hamiltonian l|2.1|) has a wide range of 
applicability, such as the Bogoliubov-de Gennes Hamil- 
tonian in superconductivity, graphite systems |l7j . fer- 
roelectricity of organic materials and perovskite oxides 
|2l| . and the slave boson mean field theory for spin liq- 
uid states, say. 

By the Fourier transformation c x = 
TV -1 / 2 X^feeBz elkX(Z k wnere the summation over k 
extends over the 1st Brillouin zone (Bz), k = 2irn/N 
(n = l,...,N), the Hamiltonian in the momen- 
tum space is given by TL — X^feeBz c I^(^) c fc' wrtn 
H ( k ) ■= J2x e ~ ikxH ( x )- If we introduce an "off-shell" 
four-vector R^= 01 ' 2 - 3 (k) e K by R°(k)TR 3 (k) := t±(k), 
— R 1 (k) + \R 2 (k) := A(fc), we can rewrite the Hamilto- 
nian in the momentum space as 

TL = 4R^k)a^c k , (2.2) 

feeBz 

where <x M = (00, — f) with 00 = 12- 

Observing that R^a^ is diagonalized by the same eigen 
vectors as those of RiOi = R- a (but with different eigen 
values), normalized eigen states v± for R^a^ are given 
by, when R is not in the Dirac string, (R 1 ,^, 2 ) ^ (0,0), 

m 

= vmm ( ±r-r 3 ) ' (2 - 3) 

where R = \R\ (should not be confused with i?°), and 
the eigen value for v± is E± = R° =F R. The Hamiltonian 



is then diagonalizcd as H — "}2 k cx\dia.g{E + , E-)kOL k , 
where c L ^ = (v r7 ) L a (Tt i c - As we assume there is a finite gap 
for the entire Brillouin zone, E+ > E-, y k g Bz. The 



vacuum |\I>) is the filled Fermi sea 111/) 



n 



JO). 



The Berry phase can be defined through the expecta- 
tion value of the twist operator: 



exp 



.2tt x 
N ^ 



(2.4) 



where n x is the electron number operator at site x, 
n x = J2 L c t;.i. c x.c- This operator twists the phase of wave 
functions along the x-dircction over the wide length scale, 
N. If we use the S z component of spin operator, say, in- 
stead of n x , we can define the twist operator in spin sys- 
tems in a similar fashion. The twist operator has been 
used to to characterize low-dimensional quantum systems 
|l6| and to describe macroscopic polarization of insula- 
tors 0, say. 

For the Fermi-sea |\1/), the expectation value of the 
twist operator is calculated as 

= {-l) N+1 exp [i 7 - $JN + 0(l/N 2 )} ,(2.5) 

where the Berry phase (Zak's phase) 7 is given by a line 
integral of the gauge field A(k) over the ID Brillouin zone 

(Bz) [HQEilll, 



\A x {k) := {v-{k)\ — |«_(fc)>, 



7 := 1 



2tt 



iA x (k) dk . 



(2.6) 



For the Fermi-sea |^) derived from the Hamiltonian 
(|2.2|l , 7 is simply equal to half of the solid angle sustained 
by a loop defined by R(k) in .R-space (Fig. Q 

On the other hand, the 0(1/N) correction to ln^l^) 
is real and given by the integral of the quantum metric 
gxx(k) over the Bz p3j . 

9xx(k) := Re{dkV-\dkV-) - (dkV-\v-){v-\dkV-), 

tic ■= * / 9x X (k)dk. (2.7) 
Jo 

The localization length £; oc plays a similar role to £ C orr 
and is known to be related to most localized Wannier 
states in an insulating phase. p3 | 



A. Truncated correlation matrix and its zero 
modes 

We next partition the system into two parts, A 
{ x x I .... . N A } and B = {x \ x = N A + 1, ... , N B } 
with Na + Nb = N, and ask, with the von Neumann 
entropy Sa, to what extent these two subsystems are en- 
tangled. Instead of directly tracing out the subsystem B 
following the definition we can make use of corre- 

as shown in Ref. 




FIG. 1: (Left) The loop denned by a three components vector 
R(k) associated with the Hamiltonian in momentum space 
[Eq. 12.H l. (Right) Loops for chiral symmetric Hamiltonians. 



|24l |. From the entire correlation matrix, we extract the 
submatrix {C L \{x^y)} x ,yeA where x and y are restricted 
in the subsystem A. The entanglement entropy is then 
given by 



•5'/ 



E 



[Ca 



lnC a + (l-Ca)hi(l-Ca) , (2.8) 



where £ a are the eigen values of the truncated correlation 
matrix {C L \(x - y)} x , y eA- 

With the whole set of the eigen values {E±(k)} and 
eigen wavefunctions {v±(k)} (Eq. (|2.o|l ) in hand, the cor- 
relation matrix C LX (x~y) = N' 1 £ fceBz e~ ik ^^C iX {k) 
is calculated exactly as 



[n M (fcK 



(2.9) 



lation matrix C L \(x — y) := 



c (z, t ) c (y,A). 



where we have introduced an "on-shell" four- vector by 
n M = (1, R/R). It should be noted that a set of Hamilto- 
nians can share the same ground state wavefunction and 
thus the same correlation matrix. 

The basic idea we will use to discuss the entanglement 
entropy is to think that the correlation matrix C(x — y) 
defines a ID "Hamiltonian" with PBC. This "Hamilto- 
nian" (let us call it the correlation matrix Hamiltonian 
or the C-Hamiltonian for simplicity) has the same set of 
eigen wave functions as the original Hamiltonian but all 
the eigen values are given by either 1 or 0. The range of 
hopping elements in the generated system is order of the 
inverse gap of the original Hamiltonian. I.e., if there is a 
finite gap, the C-Hamiltonian is local (short-ranged). 

Now, all we need to know is what energy spectrum the 
C-Hamiltonian will have when we cut it into two parts, 
defined by A and B. This is the same question asked in 
Ref. Ill in which a criterion to determine the existence of 
zero-energy edge states is presented. There are two types 
of eigen values in the energy spectrum of the truncated 
C-Hamiltonian in the thermodynamic limit Na — * oo- 
Eigen values of the first type are identical to their coun- 
terpart in the periodic (untruncated) system. On top of 
it, there appear localized boundary states whose eigen 
values are located within the bulk energy gap. Since the 
eigen values that belong to the bulk part of the spec- 
trum arc cither 1 or 0, they do not contribute to the 
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entanglement entropy as seen from Eq. (|2.8() whereas the 
boundary modes do. 

The question is then how many boundary states ap- 
pear and with what energy when the system is truncated. 
As suggested from the the KSV formula in macroscopic 
polarization, the non-vanishing Berry phase of the filled 
band of the C-Hamiltonian implies the existence of states 
localized near the boundary. Here, note that the Berry 
phase for the generated system (C-Hamiltonian) is iden- 
tical to that of the original system, since the original and 
generated Hamiltonians share the same set of eigen wave 
functions. 



B. Dimerized limit 

To know the number of localized states that appear in 
the spectrum and the energy eigen values thereof is, in 
general, a difficult task. In this subsection, we consider a 
limiting situation in which the localization length in Eq. 
I|2.7(l is small £ CO rr "C 1 and the band width of the energy 
spectrum is zero. More precisely, let us consider the case 
in which the correlation matrix is given by a four- vector 



(Eq. <E3) with 



R(k) = (-Acosfc, -Asinfc,£), 



(2.10) 



where A,£ £ R. There is a family of Hamiltonians hav- 
ing this correlation matrix which includes the following 
"dimerized" Hamiltonian, 



^ C Xi°XL + AC 1 + 1, + C X,- + h -' 



(2.11) 



The inverse Fourier transformation of Eq. I|2.10|l gives the 
correlation matrix in the tight-binding notation, 



E 



eg - 

2R 



Cr,., Crf 



A 

2R ( 



h.c. 



(2.12) 



This C-Hamiltonian can be diagonalizcd for both periodic 
and truncated boundary conditions by introducing the 
"dimcr" operators via d 



t 

±,x+) 



(4 i+ ±4 +1) _)/\/2. (See 

also Appendix.) The truncated C-Hamiltonian has (N — 
l)-fold degenerate eigen values £ = 0, 1, and two eigen 
values C = (1 ± |)/2 that correspond to edge states. 
The entanglement entropy (in the thermodynamic limit) 
is then computed as 



7 In 7 



(2tt-7), 2tt 
m ■ 



(2.13) 



2 2tt 2tt 2tt 2tt 

where the Berry phase 7 for the correlation matrix (|2.1U|) 
is 7/7r = l—£/R. In the two extreme cases, £ = and £ — > 
±00, we have S A (£ = 0) = 2 In 2 and S A (£ -> ±00) = 0, 
respectively. The entanglement entropy in the present 
case is a convex function with respect to 7 e [0, 27r] and 
the maximum is achieved when 7 = tt whereas two min- 
ima are located at 7 = 0, 2-7r. 



C. Case of 7 = n with chiral symmetry 

Although the formula l|2.13|) clearly shows the rela- 
tion between the Berry phase and the entanglement en- 
tropy in a specific limit, it is rather difficult to extend Eq. 
(|2.13(l to more generic situations. However if we impose 
a discrete symmetry implemented by a unitary particle- 
hole transformation, so-called chiral symmetry, on the 
C-Hamiltonian, it is possible to make a precise predic- 
tion for the number of boundary states that has an eigen 
value C = 1/2, following the same line of discussions in 
Ref.Gj 

When the system respects the chiral symmetry, we can 
find a unitary matrix that anti-commutes with the one- 
particle Hamiltonian. For this case, n(k) is restricted to 
lie on a plane cutting the origin in _R-space, which in turn 
means that the Berry phase for the lower band of Ti. is 
equal to rm (n £ N). (Fig. [TJ 

When n is odd, we can show that there are at least a 
pair of boundary modes at £ = 1/2, one of which is local- 
ized at the left end and the other at the right. [2f| (The 
system with 7 = ±ir is, in a sense, "dual" to that with the 
vanishing Berry phase where there is no boundary state. 
See Appendix.) Basically, this is because, when n is odd, 
it is always possible to deform the C-Hamiltonian into 
a "reference" one without closing the bulk energy gap 
and without changing the Berry phase. The reference 
C-Hamiltonian is similar to the dimerized example (|2.12|) 
for which one can exactly show the existence of n pairs 
of edge modes at £ = 1/2. In the course of deformation, 
the edge modes present in the reference C-Hamiltonian 
can move away from Q = 1/2. However, due to the chi- 
ral symmetry, the edge modes can escape from Q = 1/2 
only in a pair wise fashion, i.e., an edge state localized 
on the left /right must always be accompanied by the one 
localized on the same end and with the opposite eigen 
value with respect to Q = 1/2. When n is odd, a pair of 
edge modes ( one for each end ) cannot have its partner 
and hence we are left with at least one edge mode per 
boundary located exactly at £ = 1/2. See Ref. ^3] f° r 
more detailed discussions. 

Then, the lower bound of the entanglement entropy is 
given by 



S A > 



In- 



1 



In. 



1 



2 In 2. 



(2.14) 



This lower bound is equal to the entanglement entropy 
contained in a dimer for each end of the original model, 
which is consistent with the fact that the origin of the 
boundary states discussed above can be traced back to 
dimers in the reference Hamiltonian to which a given 
target Hamiltonian is adiabatically connected. 

There can be other contributions from boundary states 
that are not connected to a dimer in the above sense. 
Indeed, as we will explicitly demonstrate below, this kind 
of boundary modes proliferate as we approach a quantum 
critical point whose number grows as ~ A(c/6) ln£ corr /a, 
and finally gives rise to the logarithmic divergences at the 
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Dimerization 

FIG. 2: The energy spectra of (a) the Hamiltonian TC with 
open ends, (b)the truncated correlation matrix C, and (c) the 
matrix <5> (see Sec. Mil as a function of the dimerization pa- 
rameter <j) G [—1,1] f° r the SSH model. Both energy and 
dimerization are measured in the unit of the hopping ampli- 
tude, t. (d) The entanglement entropy of the SSH model. 



Under the PBC, the SSH Hamiltonian can be diagonal- 
ized as Eq. I|2.2|) with R x (k) = — 1 — — (1 — 0)cosfc, 
Ry(k) = (-1 + 0) sin k, R z {k) = 0. 

For G [—1,0), the Berry phase is given by 7 = it 
whereas for G (0, 1] 7 = 0. These two phases are 
separated by a quantum phase transition at = 0. Fol- 
lowing the discussion in Ref. 01; there is at least pair 
of boundary states for <p G [—1,0) when we terminated 
the system. Indeed, for the numerically computed energy 
spectrum of the SSH model with open ends (Fig. |2J(a) ) 
for <fi G [—1, +1], there is a pair of edge states in the bulk 
energy gap when G [—1,0). 

The entanglement entropy is calculated by diagonal- 
izing the C-Hamiltonian. The energy spectrum of the 
C-Hamiltonian with open ends is shown in Fig. 0-(b). 
Again, there is a pair of boundary states for G [—1,0) 
and for this case, Sa is bounded from below as Sa _ 
2 In 2. ( Fig. 13(d)) When we approach the transition 
point = 0, some bulk eigen values turn into the bound- 
ary eigen values and they give rise to extra contributions 
other than the zero-energy boundary states. Similar be- 
havior of the entanglement entropy is discussed for the 
quantum Ising chain in transverse magnetic field, where 
the 2 In 2 entropy originates from a Schrodinger cat state 
composed of all spin up and down configurations. 



critical point, |2| 

Note also that our discussion here does not apply 
gapless systems since the matrix elements of the C- 
Hamiltonian are long-ranged in this case. 



III. CONNECTION TO A KINK OPERATOR 

We have seen that bipartitioning the system corre- 
sponds to an introduction of a sharp "boundary" (in- 
terface). In this section, we will realize it by a non-local 
operator, a kink operator 



D. Example : the Su-Schrieffer-Heeger model 

As an example, let us look at a situation in which two 
phases with the Berry phase 7 = ir and are connected 
by a quantum phase transition point. Physically, such 
an example is provided by the Su-Schricffcr-Hcegcr (SSH) 
model for a chain of polyacetylene. The ID tight-binding 
Hamiltonian for the SSH model for a chain of polyacety- 
lene is given by H = t(-l+ (-l)Vi) {4 c i+i +h.c.) 
[26j where 0j represents dimerization at the i-th site, 
and an alternating sign of the hopping elements reflects 
dimerization between the carbon atoms in the molecule. 
Here, we treat the lattice in a classical fashion and ne- 
glected its elastic (kinetic) energy. Taking 0j = = 
const., t = 1, and defining a spinor at i = 2i — 1 



by c x = (cf,c<. 
(N = Ni/2) 



N 



the Hamiltonian can be written as 



-(1 



x \ -(1 + • 



c x+1 + h.c. 



(2.15) 



cxp 



i <p(x)n x 



r/ 1 =?/ 



where 



0, x G A, 

7T, X G B. 



(3.1) 



(3.2) 



The geometric mean of this kink operator is the twist 
operator. p7j 

The kink operator attaches a phase factor ip(x) for the 
fermion operators at site x, 



(3.3) 



Thus, if we introduce the reduced density operator 
through 



p A := - fo|*)<*tot + |*)(*|] , (3.4) 

the matrix elements tr [c XjL c y x Pa] are vanishing when- 
ever x £ A and y G B and vice versa, whereas they 
coincide with the correlation matrix C L \(x — y) when 
x,y G A. Unlike pA, the matrix elements tr [c£ jt c \Pa\ 
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are non-zero even for the B subsystem. This "padding" 
does nothing however. 

In the following, we will discuss the expectation value 
of the kink operator (*f f |?7|\I r ) with respect to a given 
ground state wave function Y$>) which is related to the ex- 
pectation value of pa as (^\pa\^) = \ [|(*|?7|*)| 2 + l]. 
As we will see the vanishing of ( \I/ 1 77 1 \I/ ) is closely related 
to a In 2 contribution to Sa discussed in the previous sec- 
tion. This can be understood intuitively as follows. Clas- 
sical wave functions can be written as a product state 
and are rather insensitive to the kink operator. Thus, 
the ground state with the kink operator inserted r]\^>) 
has a large overlap with the original ground state \^f). 
On the other hand, the kink operator destroys dimers if 
the Berry phase of the ground state is ir x (odd integer). 
As a consequence the overlap ('Jlryl'J') is very small in 
this quantum phase, which in turn suggests that quasi- 
particles that constitute the continuum spectrum above 
the ground state can be interpreted as a kink created by 
77. Thus the kink operator is capable of distinguishing the 
quantum phases with different entanglement properties. 

To put the above statement in a quantum information 
perspective, remember the reduced density matrix pA is 
in general in a mixed state: 



PA = 0)(* n (g)0|, 



(3.5) 



where \$> n ) belongs to the subsystem A, and J2 n Pn = 1- 
When the wavefunction \^f) happens to be a completely 
entanglement-free, product state, I*) = \^a) <8> \^b), 
the reduced density matrix pa is in a pure state, i.e., 
Pn^i = 0, pi = 1, l^i) = \^a)- On the other hand, 
when |\&) is highly entangled, taking partial trace over 
the B subsystem generates many pure states |\& n ) with 
non-zero weight < p n < 1. How far a given state \^f n ) 
from a product state can then be measured by taking the 
expectation value of the reduced density matrix pa- 

(*|Pa|*) = ^p„(*|*„®0)(*„®0|*). (3.6) 

n 

Clearly, it is equal to one when \^> n ) is a product state 
whereas it is expected to be less than one for entangled 
states. 

In the following subsections, we will establish that in 
an insulating phase the expectation value of the kink op- 
erator is zero in the thermodynamic limit when the Berry 
phase is 7r x (odd integer) whereas it is finite otherwise. 



A. Expectation value of the kink operator as a 
determinant 

The computation of the expectation value of the kink 



operator for a Fermi-Dirac sea |\&) = JlfceBz - /JO) S oes 
as follows. In the momentum space, the phase attach- 
ment transformation (|3.3|) reads 



fq C k- 



t t 



(3.7) 



where we introduced the Fourier components of e ltp ^ by 
= f{x) = £ /?e i<^ ( 3. 8) 

qSBz 

with q = 2im q /N (n q e N) and 



1-e" 



\ _ e -'\2-Kn q /N 

r 4 

= { 1 - e ~ i27vn i/ N ' nq 



1,3, 
0,2, 



I 0, iiq 

In a basis that diagonalizcs the Hamiltonian, 



N 
N 



(3.9) 



5>£,5 M ,(3.10) 



where a 27V x 2N matrix S 



'(fct)(fe'A) 



E 

q 



(fa)(fc'A) is S iven by 
/* [v\k-q)v(k)]^S k - q , k ,,(3.11) 



and v^(p) = (v^_(p),v <! _(p)) . 

The expectation value of the kink operator with re- 
spect to \^f) is then represented as the determinant of 
N x N matrix S, 



'(*-)(*'-)' 
(^M*) = det 



'(fc-)O'-) 



If we define the "hopping" elements i Pi9 through 



tk,k- 



— u 



(3.12) 



(3.13) 



the matrix S(k-)(k'-) hi Eq. (|3.12|) can be represented by 
a tight- binding Hamiltonian as, 



k,k' 



— ^ 2^ h t k,k-q a \ a k-q 
k q 



(3.14) 



where aj (a p ) represents a fcrmionic creation (annihila- 
tion) operator defined for p £ Bz. This Hamiltonian can 
be interpreted as describing a quantum particle hopping 
on a ID lattice. Note that the gauge field A x (k) and the 
metric g X x(k) are related to the phase and the amplitude 
of the nearest neighbour hopping elements t k fc _ 27r /jy 
respectively. The hopping matrix t k k _ q is generically 
non-local. Also, since the kink operator introduces a 
sharp boundary in the real space, the dual Hamiltonian 
is highly non-local in fc-space. 

It is evident from Eq. (|3.12|) that the vanishing of 
(^Irjl^) is equivalent to existence of zero modes in the 
spectrum of the 5-Hamiltonian. As we will see below, the 
spectrum of the 5-Hamiltonian is pretty much similar to 
that of the C-Hamiltonian: away from a critical point, the 
spectrum is gapped and all the eigen values are close to 
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cither +1 or —1, except a few eigen values in the gap that 
reflect the Berry phase if it is non-trivial. If the Berry 
phase is ir x (odd integer), there are exact zero energy 
eigen modes. When we approach a critical point, eigen 
values proliferate around zero energy. Roughly speaking, 
the entanglement entropy takes into account the distri- 
bution of all the eigen values of 5, whereas the kink 
operator only takes into account the products of all the 
eigen values. 



arg fq 


+7T/2 

^^^--^ q 


— TV 






-tt/2 


+ 7T 



B. "Chiral symmetry" and "time-reversal FIG. 3: 

symmetry" 1,3,5,- 



arg f q = arg 
• ,N-1. 



= -arg 



for n a = 



The 5-Hamiltonian has a chiral symmetry. It directly 
reflects our bipartitioning the original system and has 
nothing to do with the chiral symmetry in the original 
system. Indeed, from Eq. I|3.9[) . one can see that a k with 
k odd/even are connected to a k > with k' odd/even only. 
All the eigen states in /.'-space are connected to their 
partner with the opposite energy via 



k 



2lTTl k 

N ' 



which in turn means in the real space 



O-x+Na — a x- 



(3.15) 



(3.16) 



When the original system respects the chiral symme- 
try (not to be confused with the chiral symmetry above) , 
all the single particle wave function ip(k) of iS in fc-space 
can be taken to be real by a suitable rotation in U-space. 
[However when the Berry phase is 7 = it x integer, this 
comes with a price to have a Dirac string that intersect 
R(k).] The ability of taking all "hopping" elements tk,k' 
to be real induces an additional "time-reversal symme- 
try" to 5-Hamiltonian; the phase associated with f q can 
be removed by a simple gauge transformation, 



CLk 



3 +ifc/2-ifciV A /2 



Ofc. 



(3.17) 



[See Fig. [3J] Thus, we can take all the matrix elements 
f q t k k _ q in the 5-Hamiltonian to be real. Furthermore, 
when we go back to the real space, this "time-reversal" 
invariance implies a parity symmetry with respect to an 
inversion center xo = —Na/2 + 1/2. To see this, we first 
note that all the one-particle eigen states of S can be 
taken real in the basis {6^,6 p }; the 5-Hamiltonian can 
be diagonalized as S = J2 n en ^n^n with 



(3.18) 



where <f> n (p) is a eigen wavefunction which is real. Since 
the basis {al.,a x } and {<^,c£ n } are related through 

a - = E4fE eifc( ^ 1/2+AU/2) ^( fc H-( 3 - 19 ) 



the real space eigen wavefunctions ip n (x) in the basis 
{ a x> a x) are given by 



1 e ik(x-l/2+N A /2) 



Mk), (3-20) 



from which one can see ip n (x) satisfies 

[il> n (x)] m = ^ n (-x + l-N A ). 



(3.21) 



I.e., the wave function amplitude is parity symmetric 
with respect to xo = —Na/2 + 1/2. 

This time-reversal symmetry, which is plays an impor- 
tant role for the vanishing of the expectation value of the 
kink operator. Indeed, it is this symmetry which guar- 
antees existence of zero-modes of S. 



C. Existence of zero- modes 

The argument that tells us the existence of zero modes 
for the 5-Hamiltonian is somewhat similar to the "proof" 
of the existence of zero modes for the C-Hamiltonian in 
that we consider a adiabatic change of the Hamiltonian. 
The major difference comes from the fact that the chiral 
symmetry in the 5-Hamiltonian is implemented as a kind 
of time-reversal symmetry as we discussed before. 

We first establish that there is a pair of zero modes for 
S when we take |^) as the ground state of the dimerized 
Hamiltonian l|2.10|l with the chiral symmetry. The hop- 
ping elements in S arc computed from the overlap of the 
Bloch wave functions as 



(v±(p)\v±(q)} 
1 



A 2 e i(p_?) + R 2 T 2i?£ + £ 5 



2R(RtR 3 ) 
The 5-Hamiltonian is then diagonalized as 
1 



(3.22) 



S 



2R(R-R 3 ) 

x £ [A 2 f(x + 1) + (R- €) 2 /(x)] 4a x . (3.23) 
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We see that there are two mid-gap states with energies 
±£/R. Especially when £ = 0, there are a pair of zero 
energy states localized at the interfaces. 

We then change the Hamiltonian in a continuous fash- 
ion in such a way that (i) it respects the chiral symmetry 
during the deformation, and (ii) it does not cross the 
gap closing point (the origin of i2-space). During this 
deformation, the Berry phase of the ground state wave- 
function is always kept to be %. As already discussed, 
we can take all the Bloch wave functions to be real and 
there is a "time-reversal" symmetry. 

One can see that the zero modes never escape from E = 
as it constrained by the time-reversal symmetry, which 
is nothing but the parity invariance with respect to xq = 
—Na/2 + 1/2. First note that since the 5-Hamiltonian 
in fc-space is non-local, it is short-ranged (quasi-diagonal) 
in the real space. Thus, if we take the thermodynamic 
limit N — ► oo, states that appear between the gap are 
spatially localized near the interfaces located x = 1/2 
and x = Na + 1/2, which separate the system into the 
two subsystems. 

During the deformation, the two localized states, which 
located at x = 1/2 and x = Na + 1/2, respectively, can 
in principle go away from E = 0. Due to the "chiral 
symmetry" of the 5-Hamiltonian, if one goes up from 
E = 0, the other must be go down. However, if there 
is the "time-reversal symmetry", each eigen state must 
be invariant under the space inversion with respect to 
—Na/2 + 1/2. In order for the localized states to satisfy 
these two conditions, both of them must be located at 
E = 0. 

As an example, the spectrum of the 5-Hamiltonian for 
the SSH model is presented in Fig. [3(c). The spectrum 
is almost identical to that of the C-Hamiltonian and a 
pair of zero modes persists for the entire quantum phase 

0e [-1,0). 

IV. 2D SYSTEMS WITH THE 
NON- VANISHING CHERN NUMBER 

As far as we consider translational invariant systems, 
the above ID discussions still apply to higher dimensions. 
When a d-dimensional translational invariant system is 
bipartitioned by a (d— l)-dimensional hyperplane, we can 
perform the (d — l)-dimcnsional Fourier transformation 
along the interface. The Hamiltonian is block-diagonal 
in terms of the wave number along the interface fe|| , TL —: 
Sfe ^(fyl ) j where TC(k^) is a ID Hamiltonian for each 
fc|i-subspace. Then, the previous discussion applies to 
each H(feii). As an example of a 2D two-band system, 
let us consider 2D chiral p-wave superconductor (p-wave 
SC) defined by 

H = 4 ( _a -f ) Cr+£ + h,c ' 

+c - ( il *-t ) Cr+ y + Lc - + c - ( o ) Cr ' (4 - x) 



where the integral index r runs over the 2D square lattice, 
x = (1, 0), y = (0, 1), and t, A, \x £ R. For simplicity, we 
set t = A = 1 in the following. The chiral p-wave SC has 
been discussed in the context of super conductivity in 
a ruthenate and paired states in the fractional quantum 
Hall effect. M, HJ El El E3 There are four phases 
separated by three quantum critical points at \x = 0, ±4, 
which are labeled by the Chern number Ch as Ch = 
> 4), Ch = -1 (-4 < n < 0), and Ch = +1 (0 < 
H < +4). The non-zero Chern number implies the IQHE 
in the spin transport. |3oT | 

The energy spectrum of the family of Hamiltonians 
Ti.(k y ) parametrized by the wave number in y-direction, 
k y , is given in Fig. 0J-(a),(b). There are branches of edge 
states that connects the upper and lower band for phases 
with Ch = ±1. These edge states contributes to the 
entanglement entropy. The energy spectrum of the C- 
Hamiltonian with open ends is shown in Fig. 0J-(c),(d). 
The corresponding entanglement entropy is also found 
in Fig. 01(d) for several values of the aspect ratio r = 
N y /N x . We can see that for small r, the entanglement 
entropy shows a cusp-like behavior at quantum phase 
transitions whereas for larger value of r, the cusp is less 
eminent. 

This behavior can be understood as a dimensional cross 
over of the scaling behavior of the entanglement entropy 
between ID and 2D. For small r, the entropy behaves 
ID-like and the cusp is a reminiscence of the logarith- 
mic divergent behavior Sa ~ In Na of the pure ID case. 
Q On the other hand, for r close to unity, the entropy 
exhibits a 2D behavior. In the pure 2D limit (r = 1), not- 
ing that the band structure at the critical points [i = ±4 
consists of one gapless Dirac fermion, the entanglement 
entropy scales as Sa = aN y —(3N y /NA where a, (3 is some 
constant. (See Appendix IbI) Notice that unlike the case 
of a finite Fermi surface |3J;lll|i Sa/N v is constant for a 
Dirac fermion. 

An interesting and direct application of the present 
section is the entanglement entropy of 2D d-wave super- 
conductors and carbon nanotubes. In these systems, dif- 
ferent ways of bipartitioning the system lead to different 
amounts of the entanglement entropy. ^} 



V. CONCLUSION 

In this paper, we have identified two types of con- 
tributions to the entanglement entropy, i.e., one from 
the boundaries of the system created by taking the par- 
tial trace and the other from the bulk energy spectrum. 
The contribution from the boundaries is controlled by 
the Berry phase and hence we can make use of some 
known facts on the "bulk-boundary correspondence" to 
compute the entanglement entropy. Especially, we have 
obtained the lower bound of the entanglement entropy for 
ID systems with discrete particle-hole symmetries. Intu- 
itively, this means that when the Berry phase is zero, 
the ground state wave function is very close to a sim- 
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FIG. 4: The energy spectrum (measured in the unit of the 
hopping t — 1) v.s. k y £ [0, 2tt) for the 2D p-wave SC with 
boundaries. The chemical potential is n — —5 (a) and —3 
(b) and t = A = 1. The corresponding spectra of the C- 
Hamiltonian are shown in (c) (/i = —5) and (d) (^ = —3). 
The entanglement entropy of the 2D chiral p-wave SC as a 
function of fj, is presented in (e).n The aspect ratio r — N y /N x 
is r = 1/2, 1/3, 1/4, 1/8 from the bottom at fi = -4. 



pie product state, and there is not much entanglement. 
Thus, ground states with non-trivial Berry phase can be 
said to be more entangled in general. 

Recently, it has been revealed that the Berry phase 
manifests itself in the semiclassical equation of motion 
|36j . the density of states [53, and the anomalous Hall 
effect, etc. One can put the Berry phase correction to 
the entanglement entropy in the catalog. 

One of the main massages of this paper is the superi- 
ority of the entanglement entropy to conventional corre- 
lation functions of local operators to describe quantum 
phases. Indeed, we clarified that the entanglement en- 
tropy is related to non-local operators ; the twist operator 
and kink operator. The bulk contribution to the entan- 
glement entropy is related to the localization length (cor- 
relation length) which is the real part of the logarithm 
of the expectation value of the twist operator and can 
be expressed by the quantum metric [231 ]. On the other 
hand, the edge contribution is tied with the imaginary 



part and to the Berry phase. [See Eqs. i|2.4[) to i|2.7|) .] 

We have also made a connection between the entangle- 
ment entropy and the kink operator. It is known that 
several phases of ID strong correlated systems (such as 
the Haldane phase) can be described by these non-local 
operators. Another connection of the entanglement en- 
tropy to a some sort of non-local operator can be also 
seen in a recent proposal of a holographic derivation of 
the entanglement entropy. [38l | 

Thus, the entanglement entropy can potentially be 
very useful to detect several quantum phases that need a 
more subtle way of characterization than classically or- 
dered phases. For example, the entanglement entropy 
can be applied to several types of spin liquid ground 
states which are speculated to be described by some kind 
of gauge theories. Indeed, for gapped phases of topolog- 
ical orders, this direction has already been explored to 
some extent [3^, . 

However, in order to push this direction further, we 
still need to deepen our understanding of the entangle- 
ment entropy. For example, extensions to multi-band 
systems, especially to the case of completely degenerate 
bands might be also interesting in which we need to use 
the non-Abelian Berry phase to characterize the system. 
[4l| It is also interesting to investigate if the Berry phase 
of quantum ground states can be captured by other t ype s 
of entanglement measures such as the concurrence |42| . 
Finally, among many other questions, we need to con- 
sider how we can measure the entanglement entropy in a 
direct fashion. l43l 
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APPENDIX A: THE DUAL BERRY PHASE 

In this appendix, we introduce the dual Berry phase 7 
in the ID two-band Hamiltonian Eq. I|2.1(l . If we impose 
the chiral symmetry, in a quantum phase with 7 = 
the dual Berry phase is given by 7 = — it whereas when 
7 = — 7T, 7 = 0. Thus a quantum phase is characterized 
by both 7 and 7. 

It is in spirit similar to the dual order parameter (dis- 
order parameter) in the quantum Ising spin chain. The 
ID quantum Ising model in a transverse field has two 
phases : ordered and disordered phases. It is known that 
the entanglement entropy is Sa > 2 In 2 for the former 
and Sa > for the latter. On the other hand, they 
are related to each other by the Kramers- Wanicr duality 
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and hence one may argue that they arc essentially equiv- 
alent. Why is the entanglement entropy in the ordered 
phase is larger than that in the disordered phase ? The 
reason is that Kramers- Wanier duality transformation is 
a non-local transformation and does not leave original 
bipartitioning invariant. 

Similarly, the duality that we will introduce momen- 
tarily connects two different Hamiltonians with different 
Berry phase and entanglement entropy ft is possible 
since it is a transformation that changes the way of la- 
beling of sites and hence bipartitioning. 

Let us first introduce dimer operators by 



±,x + 



1 

3 



,x+l 



(Al) 



When written in terms of the dimer operators, the Hamil- 
tonian (|2.1|l reads H = \ J^x' ^ _l i H x - X <d i , where 

the new hopping matrix elements H x - X ' are some func- 
tion of the original ones H x ^ x >. For simplicity, we focus 
on the case of particle-hole symmetric (t + = — f_) and 
translational invariant systems. In the momentum space, 
the Hamiltonian is given by Ti 
with a 3D vector R(k) given by 



EfcGBz d k R ( k ) ■ ad k 



Rx{k) = R z (k), 

Ry(k) = sm(k)R x (k) + cos(k)R y (k), 

R z (k) = cos(k)R x (k)-sin(k)R y (k). (A2) 

We define the dual Berry phase 7 as the Berry phase for 
the dual 3D vector R(k), 



7 - 



dk(v(k)\ — \v(k)) 



dk- 



XY - YX 
V 2R(R - Z) ' 



(A3) 



where X = dR x (k) /dk, etc. Rotating R around R y - 

axis as (R x ,R y ,R z ) — > (R z ,R y , —R x ), and noting Y X — 

XY = XY - YX + X 2 + Y 2 , the dual Berry phase is 
thus given by 



7 = —7 — 7T 



dk — 

2R 



(A4) 



Especially if we impose the chiral symmetry, R is re- 
stricted to ATY-plane and thus 



7 



-7 — 7T. 



(A5) 



APPENDIX B: ENTANGLEMENT ENTROPY 
FOR A DIRAC FERMION IN 2D 

In this Appendix, we estimate how much entangle- 
ment entropy is carried by a gaplcss Dirac fcrmion. It 
is an interesting question since a Dirac fermion is just 
in-between of a fully gapped system and a system with 
a finite Fermi surface; for the former the entanglement 
entropy satisfies the area law, whereas for the latter there 
is a log-correction to the area law. 0, [3f| 

Unlike the case of a finite Fermi surface, the entan- 
glement entropy divided by N y , Sa/N v is constant for 
a system with Dirac fermions as can be seen as follows. 
The energy spectrum close to a gap closing point in the 
Bz, fc' ', is linear and so is the mass gap as a function of 
k y , m(k y ) oc ky — k y °\ The known result for a massive 
ID system tells us each k y contributes to the entangle- 
ment entropy by ~ \nm(k y ) . If the length Na of the 
subsystem A in cc-direction is finite, we expect the contri- 
butions from those k y with m{k y )^ 1 > Na are given by 
~ IjiNa, instead. Then, the entanglement entropy can 
be evaluated by summing over the entanglement entropy 
for each ID system with a fixed k y , 

4= E fw*.r l + E Ik*.*. 

fc B >0 fc y >0 

(Bl) 

Converting the summation to the integral, we see that the 
entanglement entropy behaves as Sa = aN y — (3N v /Na 
for a single Dirac fermion where a, /3 is some constant. 
Hence Sa/N v is finite. This crude approximation is actu- 
ally overestimating the entropy, but it is enough to derive 
essential features of the entropy. For more detailed anal- 
ysis using the entropic c- function, see Ref . 01 ■ 
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